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It is a central question in quantum thermodynamics to determine how much work can be gained 
by a process that transforms an initial state p to a final state a. For example, we might ask how 
much work can be obtained by thermalizing p to a thermal state a at temperature T of an ambient 
heat bath. Here, we show that, for different sets of resource-theoretic thermodynamic operations 
involving catalysis, the amount of work is characterized by how reversible the process is. More 
specifically, the amount of work to be gained depends on how well we can return the state a to its 
original form p without investing any work. We proceed to exhibit an explicit reversal operation in 
terms of the Petz recovery channel coming from quantum information theory. Our result establishes 
a quantitative link between the reversibility of thermodynamical processes and the corresponding 
work gain. 


I. INTRODUCTION 


Quantum thermodynamics is experiencing a renaissance in which ideas from quantum information theory enable 
us to understand thermodynamics for even the smallest quantum systems. Our inability to apply statistical methods 
to a small number of particles and the presence of quantum coherences make this a challenging undertaking. Yet, we 
are now indeed able to construct very small quantum devices allowing us to probe such regimes experimentally [1—3] . 
Theoretical results studying the efficiency of small thermal machines [4-8], catalysis [9-11], work extraction [12-19], 
and the second laws of quantum thermodynamics [20, 21] have furthermore led to the satisfying conclusion that the 
usual laws of thermodynamics as we know them can be derived from the laws of quantum mechanics in an appropriate 
limit. 

Here we will be concerned with the fundamental problem of how much work is gained (or invested) by the transfor¬ 
mation of a state ps to a state as of some system S in the presence of a thermal bath. The second laws [20] provide 
general bounds on the amount of work, which are tight if ps is diagonal in the energy eigenbasis of the system. Special 
instances of this problem have drawn particular attention, such as gaining the maximum amount of work from ps by 
thermalizing it to the temperature of the surrounding bath [12], extracting work from correlations among different 
subsystems when ps is a multipartite state (see, e.g. [22]), as well as the case when as results from a measurement on 
ps [23-25]. When thinking about investing work, one of the most well studied instances is Landauer’s erasure [26], 
which is concerned with the amount of energy necessary to take an arbitrary state ps to a pure state as- 

We adopt the resource theory approach of [ .2, 27, 28], which has the appealing feature of explictly accounting for all 
energy flows. Let us first establish what we mean by gaining work. In the macroscopic world, work is often illustrated 
as raising a weight by a certain amount. An analogue in the quantum regime corresponds to raising a system from 
a state of low energy to an excited state. We will call a system that is used to store energy a battery. The simplest 
model of such a battery is given by the ‘wit’ [12], a two-level system that is normalized such that the ground state 
|0) has energy ‘O’, and the excited state |1) has energy W, corresponding to the amount of work W we wish to store. 
That is, the Hamiltonian of the battery is given by Hw = W|l)(l|. Note that raising the system from the ground state 
to the excited state not only changes its energy by W, but this energy transfer takes place in a fully ordered form: 
we know the final state and can thus, in principle, transfer all of the energy onto a third system. This is analogous to 
the macroscopic notion of lifting a weight on a string and thus does not include any contributions from heat. Clearly, 
in a specific physical system, it can be difficult to realize an energy gap of precisely the amount of work W that we 
wish to extract. However, one can imagine this two-level system as being part of a quasi-continuum battery from 
which we pick two levels with a suitable energy gap [29]. We refer to [29] for a discussion on (approximate) work in 
the microregime. When gaining work W by transforming ps to as, we are thus implementing the process 

PS ® 10X01 w -t a s ® |1)(1|vf • (1) 

The situation of investing work can be described similarly, except that instead of raising the system to the excited 
state, we draw energy from the battery by lowering it to the ground state. When investing work W in order to enable 
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a transformation, we are thus implementing the process 

Ps ® |l)(l|w —► ffs ® |0)(0|vy ■ (2) 

Within the resource theory approach to quantum thermodynamics, there are various different classes of thermo¬ 
dynamic operations. Most of the classes we will study in this paper have in common that, given a particular fixed 
temperature T, we may access a bath described by a Hamiltonian Hb and thermal state tb = sxp(—/3Hb)/Zb, where 
/3 = l/(kT) is the inverse temperature, k is the Boltzmann constant, and Zb is the partition function. To help us, we 
may also make use of catalysts in states r/Ci with Hamiltonian He = ®" =1 Hc t - Let Hg be the Hamiltonian associated 
with the system S and let U denote a unitary that acts on the system 5, the battery W, the catalysts C = ®” =1 C' i , 
and the bath B. The only unitary transformations U that are allowed are those which conserve total energy. That 
is, the allowed unitaries are such that [ U , H\ = 0, where H = Hg + Hw + He + Hb is the total Hamiltonian. The 
transformation T performing the mapping T(pg 0 |0)(0|w) = erg <8> |l)(l|vv then takes the following form 

T{wgw) = Tr cb[U(u}sw ® VCi ■ ■ ■ Vc n (3) 

for some input state cogw of the system and the battery. The resource theoretic operations we study in the paper, all 
differ in the nature of the physical restrictions imposed on the catalysts. For example, in the class of so-called catalytic 
thermal operations [20] we only allow for a single catalyst, and, in the regime of exact catalysis, we furthermore require 
that the catalyst is conserved exactly. That is, in such a case the output state on SW and the output state ty™ 1 
on C satisfies 

Tr s [U (tjj sw 0 ric ® 0 rfe t (4) 

with pc = ty™*. Whenever the catalyst is absent, the allowed transformations are a subset that is simply called 
thermal operations [12]. The other classes of resource thermodynamic operations will be discussed in Sections V 
and VI. Given that U conserves total energy, it is clear that this framework accounts for all energy flows, making it 
particularly appealing for studying quantum thermodynamics. 

How much work could we gain by transforming pg to ag using a bath of temperature T? It is clear from the above, 
that this question can be answered by asking about the largest value of W ga i n (ps —> ag) = W that can be achieved 
by a thermodynamical operation belonging to the particular class in question, e.g. catalytic thermal operations, in 
the transition of Eq. (1). The standard second law tells us that this transformation is possible only if 

F(pg®\0)(0\w)>F(a s ®\l)(l\ w ) , (5) 

where F(uigw) = Tr [Hgw^sw] ~ kT S(uigw) is the Helmholtz free energy with Hg\y = Hg + H\y and S(ujsw) = 
— Tr[wgty log uigw] 1 - When pg is diagonal in the energy eigenbasis, then the standard second law does in fact give 
necessary and sufficient conditions for the desired transformation to exist whenever the systems are extremely large 
[28] or if we allow for a slightly inexact catalysis [20]. Specifically, if an arbitrary catalyst can be used, what we 
mean by this is that the “error per particle” in the output catalyst is bounded as ||? 7 c — ?7c Ut ||i < e/log dc, where 
dc is the dimension of the catalyst and e > 0 is some tolerance [20]. Similarly, if inexact catalysis takes on the form 
of allowing small correlations in the output catalyst, only the standard free energy is relevant [10]. This is the case 
of correlating-catalytic thermal operations, which we will study in Section V. While this regime of inexact catalysis 
is—due to an accumulation of errors- highly undesirable when analyzing cyclic or quasi-static processes as is the case 
in thermal machines [29] , it is an otherwise very well motivated regime in physical implementations where an ever so 
slight error is essentially unavoidable. In the present Section this regime will be our primary focus. 

Using the fact that Tr[Uw’|0)(0|w] = 0 and Tr[.HV|l)(l|w] = W ga - ln (pg —> ag), it is easy to use (5) to obtain the 
following upper bound on the amount of work that we can hope to obtain 

W gain (p s ->■ as) < F(pg) - F(a s ) . (6) 

If pg is diagonal in the energy eigenbasis, then (6) is tight. It will be convenient to note [30] that the free energy can 
also be expressed in terms of the quantum relative entropy 

D(ps\\rs) ■= Tr[ps log ps] - Tr[ps log rg], (7) 

defined as above when supp(ps) C supp(rs) and equal to +00 otherwise. Note that here and throughout, we consider 
that the operator logarithm is evaluated only on the support of its argument. Specifically, F(pg) = kT[D{pg\\Ts) — 


1 All logarithms in this paper are base e. 
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logics], where rg = exp(— f3Hs)/Zs is the thermal state of the system at the temperature T of the ambient bath. 
Since we do not change the Hamiltonian of the system, we can hence express the amount of work in regimes where 
the standard free energy is relevant, as 


Wg ai „(ps -»■ a s ) = kTA , (8) 

where the difference A of relative entropies will play a very special role 

A = D(ps\\t s ) - D(<js\\ts) ■ (9) 

Note that when A > 0, the standard second law allows for the transformation ps —► erg. 

How about the case in which we need to invest work? In regimes where only the standard free energy is relevant to 
dictate the transition (2), a similar calculation yields that W[ nv (ps —> as) = kT(D(as\\Ts) — D(ps\\rs))- This means 
that in such regimes Wi nv (ps — > as) = — W ga i n {Ps &s) = Wgain {as —I Ps), he., the amount of energy that we need 
to invest to transform ps to as is precisely equal to the amount of work that we can gain by transforming as back 
to ps- In the standard free energy regime, we thus see that we do not need to treat the amount of work invested as 
a separate case, but rather it can be understood fully in terms of the transformation of as back to ps in which work 
can be gained. It is useful to note that for systems S that are truly small [12], or when we are interested in the case 
of exact catalysis, then this is not the case. In these situations, the standard second law needs to be augmented with 
more refined conditions [20] that lead to differences. With some abuse of terminology, we will refer to this as the nano 
regime. In place of just one free energy, the nano regime requires that a family of free energies F a satisfies 

F a (ps) > F a (a s ) , (10) 

for all a > 0. These generalized free energies can be expressed in terms of the a-Renyi divergences as 

F a (ps) = kT[D a (ps\\T S ) - log Z S \ , (11) 

where the general definition of D a 2 takes on a simplified form if ps is diagonal in the energy eigenbasis. More 
precisely, 

D a (ps\\rs) = Pj T j~ a ’ ( 12 ) 

3 

where p 3 and Tj are the eigenvalues of ps and ts respectively. The standard free energy is a member of this family 
for a —> 1. A short calculation [20] yields that in this regime 

W^ips -+ ^s) < inf kT (D a (p s \\r s ) - D a (as\\r s )) , (13) 

a>0 

W inT°(Ps -t as) > sup kT {D a {a s \\T S ) - D a (p s \\T S )) 

a> 0 

> kT (D(a s \\T S ) - D(ps\\t s )) , (14) 

where (the first) inequalities are again attained if ps is diagonal in the energy eigenbasis. In the nano-regime, it is 
thus possible that ± W-™ no . 


II. RESULT 

Gaining work. Our main result is to prove that the amount of work Wg a i n (ps —► as) > 0 gained when transforming 
ps to as with a particular class of thermodynamic operations, can be characterized by how well we can recover the 
state ps from as using a thermodynamic operation of the same class which requires no work at all. This result applies 
to all situations in which W is characterized by the standard free energy (see above). It forms a powerful link between 
the reversibility of the initial process and the amount of work drawn from it. Loosely speaking, we will see that if 
little work is obtained when transforming ps to as with a particular class of thermodynamic operation, then there 


For arbitrary states pg, we have for 0 < a < 1/2 that fA(Csll r .S') = 1 ( logTr/p/r), “] [31] and for a > 1/2, D,> (p,sj| T .s) — 

log [Tr (4 1 - q)/2q Ps 4 1 - q)/2 “)“] [32, 33], 
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exists a thermodynamic operation of the same class that can recover ps from as quite well. Or stated differently, if 
this thermodynamic operation performs badly at recovering ps, then the amount of work can be large. More precisely 

Wgain(ps -t Vs) > -kTlogF(ps,Kp^ a (a s )) , (15) 

where F(p,a) = Tr[^/ \fpa^/p] 2 is the fidelity, and 7 Z p ^. a is a reversal operation using a bath at temperature T , and 
(potentially) catalysts that takes as close to the original state ps- Wgain(ps —► Vs) is thus related to the reversal 
operation 

w gain 

Ps &s ■ (16) 

n 

This reversal operation 7 Z p ^ a is generated by using the same bath (and possibly catalysts too) as in the original, but 
reversing the global dynamics. That is, if the original operation is generated by a global unitary evolution U(-)U*, 
the reversal operation is generated by U' (-)U. Hence, no work is required to perform it. Below we also show that the 
reversal operation coincides with the Petz recovery map [34-37] from quantum information theory. 

Investing work. As outlined above, in regimes where only the standard free energy is relevant, the amount of 
work W- Inv (ps —> as) > 0 we need to invest to transform ps to as can be understood fully in terms of the work 
W g ain(o's —» ps) gained by the inverse process. This is not true in the nano regime. However, it is a nice application 

of our analysis to show that the reversal operation from the transformation as —► ps that allows us to gain work, 

nevertheless allows us to make statements about WV™ 0 . Specifically, 

WZ n °(Ps <t s ) > -kT\ogF(a s ,n^ P (ps)) , (17) 

where lZ a -> p is the reversal operation constructed from the process that transforms as to ps —from which work might 
be gained—while depositing such energy into the bath. That is, the reversal operation obtained from 

w gain 

Ps *=> &s (18) 

n 

can be used to transform ps to as at no energy cost to the battery, and allows us to make statements about W- mv : 

w inv 

Ps =1 &s ■ (19) 

n 

Needless to say, since for example the erasure of a thermal state ps = ts to a pure state as costs a significant amount 
of work, the operation TZ a ^ p cannot always achieve a very high fidelity. This is indeed, precisely what we see here 
since if the fidelity is very small, then the amount of work that we need to invest is large. We illustrate this more subtle 
application of our result in Section IV by means of a simple example of a harmonic oscillator bath. We emphasize 
that the bound on is valid for all the classes thermodynamic operations we study in this paper. 


III. PROOF FOR THERMAL OPERATIONS (NO CATALYSTS) 

We now give details of our main result for the simpler case of thermal operations (TO) without catalysts. Section V 
contains details of other, more general sets of operations. 

Let us first suppose that we can draw a positive amount of work by transforming ps to as, so that A > 0. Note 
that in regimes dictated by the standard free energy, A > 0 implies that there exists a different thermal operation 
taking ps to as without drawing any work at all [20, 21]— in this case the additional energy can be deposited into 
the bath. Let V be the energy-conserving unitary that realizes this latter thermal operation, and let (t b , Hb) be the 
thermal state and Hamiltonian of the bath, such that as = Tr B [V(ps <8> f B )V^]. Note that V acts on systems S and 
B and [ V , FIs + H B \ =0. We have the following theorem: 

Theorem 1. Let T be a thermal operation given by 

T(-)s = Tr B [V((-)s®T B )V% (20) 

where V and t b are defined above. Then it obeys the inequality 

D(ps\\rs) - £>(vs||vs) > D(p s \\1Z(a s )) > - log (F(p s ,H(a s ))), (21) 

where 7Z(-) is a recovery map, which is another thermal operation given by 

7^(•) =Tr b[V\(-) s ®t b )V}. 


( 22 ) 
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Proof. Our proof is divided into two main steps. 

Step 1: Rewriting the relative entropy difference. 

Our first step will be to rewrite A = D(pg\\Tg) — D(ag\\Tg) as an equality involving the operation V. Observe that 

D(ps\\ts) = D{ps ® t b \\t s ® t b ) (23) 

= D(V(ps®tb)V^\\V{ts®t b )V^) (24) 

= D(V(p s 0 r B )l /t ||Ts 0 t b ), (25) 

where we have the used the facts that the relative entropy is invariant with respect to tensoring an ancilla state or 
applying a unitary, and V is an energy-conserving unitary so that V(rg 0 tb)V^ = Tg 0 t b . 

For density operators pcd and 9cd such that suppler*) C supp(0c B ), it is possible to write 

D{t}cd\\9cd) — D{j]d\\0d) = T^{vcd [log tjcd — log Ocd — log Ic ®Vd + logic ® Od]) ■ (26) 

Using these two facts, we can rewrite A as follows: 

D(pg\\Tg) - £>(crs||rs) = Tr (V(ps 0 r B )U t [log V(pg 0 f|,)U t - logr s 0 f B - logo’s ® Ib + log Tg 0 I B ]) ■ (27) 

We can simplify the operator consisting of the last three terms on the right above as 

- log Tg 0 t b log ag 0 Ib + log Tg 0 I B = - log Is 0 tb - log ag 0 I B (28) 

= - log a s ®t b , (29) 

and thus conclude that 

D(ps\\t s ) - D(a 3 \\T S ) = D(y{pg®TB)V^\\as®T B ). (30) 

Hence we have that the right-hand side is equal to 

D(V(p s 0 T B )V^\\a s 0 t b ) = D(p s 0 T B \\V\crs 0 t b )V). (31) 

Putting everything together, we see that 

D{ps\\rs) - U>(o-s||ts) = D(p s 0 r B ||U t (o’s 0 t b )V). (32) 

Thus, the quantity A related to the work gain in ( 8 ) is exactly equal to the “relative entropy distance” between 
the original state ps 0 r B and the state resulting from the following thermal operation: 

a s ->■ Wl (a s 0 t b )V, (33) 

which consists of adjoining ag with a thermal state r B and performing the inverse of the unitary V. Note that this 
statement is non-trivial, since cts®tb 7 ^ V(ag 0 f B )V’t. The forward operation V can create correlations between the 
system and the bath, whereas V "1 is applied to a fresh and entirely uncorrelated bath, making it a thermal operation. 

Step 2: A lower bound using the recovery map. Due to the fact that the quantum relative entropy can never 
increase under the action of a partial trace [38, 39], we can conclude from (32) that the following inequality holds 

d (ps\\ts) - -D (os I Its) > D(pg\\TZ p ^ a (a s )) , (34) 

where 


n p ^ a (a s ) = Tr B [yl (as 0f B )U]) . (35) 

Note that this operation is a thermal operation, and requires no work. Our claim (15) now follows from the relation 
between the a-Renyi divergences [32] 

D{p\\cr) > Di(p\\a) = - logF(p, a) . (36) 

This concludes the proof. □ 


Investing work: Our result for the case where we invest rather than gain work follows from the very same analysis 
as above by exchanging the roles of ps and as: if -D(o\s||ts) — D(ps\\ts) > 0, then we could gain work by transforming 
ag and pg. The above argument yields a recovery map TZ a ^. p that lower bounds the relative entropy difference. Our 
claim of (17) now follows from (14). 
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A. Remark: Petz recovery map 


We remark that TZ is actually a special quantum map, called the Petz recovery map [34-37]. For a general quantum 
map A f and a given density operator 9, this recovery map is defined as 


a 7 (.) = e^Mey^^ATie)- 1 ^ 1 / 2 , 


(37) 


where Af' is the adjoint of the channel A7 [40]. As a consequence, we can conclude that the main conjecture from [41] 
holds for the special case of thermal operations. We show this in the following lemma. 

Lemma 1. The map TZ(-) in (35) is the Petz recovery map of the original thermal operation, provided we choose the 
9 in (37) to be the thermal state ts- 

Proof. Consider that for two density operators p and 9 and a quantum channel Af , we can associate the relative 
entropy difference D(r]\\0) — D(Af(ij)\\AF(0)) and the Petz recovery channel: 

A/'(•) = 9 1/2 ^[N(9)- l ' 2 {-)N{9)- 1/2 ]9 1/2 . (38) 


For our case, we have that 

9 = PS , 0 = t s , Af(-) = Tr s [F((-)s ® t b )V ] ], 

which implies that Af{6) = ts- Using the definition of the adjoint, one can show that 


A7t(.)=Tr s [(.) s , 


' Ib}Vt 


1/2 

B 


which implies for our case that the Petz recovery channel takes the following form: 


ti(-) = t s /2ti 'b 


1/2 (-)sr s 1/2 0 Ib]Vt 


r-V2 


- 1/2 


. 1/2 


A/ 2 


We can rewrite this as follows: 


Tr F 


(r| 0fJ)U t ['r s 2 (-) S t s 3 ® /b]V(t| 0t 2 ) 


= Tiq 


= T 17 


(rs^TB^V^Tg 2 (-)s T S 

V\t S ® T B )^T S 2 {-)sT s 2 


= Wfl[yt((.) s ®f B )y], 


5 I b ]V{t s 0 r B ) 2 
1 Ib\(ts 0 t b )^V 


where we have used that [V, ts 0 t b ] = 0. 


(39) 

(40) 

(41) 

(42) 

(43) 

(44) 

□ 


IV. EXAMPLE FOR THERMAL OPERATIONS 

Let us illustrate the reversal operation lZ a ^ p by means of a simple example. Let S' be a two-level system, with 
Hamiltonian Hs = 11>(11- Suppose that S is in the mixed state ps = Po|0)(0|s+pi|l)(l|s withp 0 € [1— e~ l3Es , l],and 

we invest work WPff 1 °(ps —> crs ) to bring the system to the ground state as = |0)(0|s. That is, we are performing 
Landauer erasure. Recall that the reversal operation associated with the lower bound for the work invested (17) 
is determined by the operation that takes as = |0}(0|s to ps without drawing any work, but instead dumping the 
resulting energy into the bath. 

For our simple example, consider a bath comprised of a harmonic oscillator H B = E n \n)(/n\ B where E n = 

nhu> 3 . Note that, for each n, the gap between n and n + 1 is constant: G = E n+ 1 — E n = two. To illustrate, let 
us consider the energy gap of the system to be equal to Es = two —an example in which Es is a multiple of two is 
analogous. 


3 We could have also written E n = (2 n + 1) (a;, which is the same after re-normalizing. For notational convenience we have subtracted 
the constant ^ uo. 













1. Transforming as to ps 


Our first goal is to find the explicit operation that takes as to ps, that has the effect of mixing the ground state 
of the system. We will see that no catalyst is required for this transformation, that is V = E7 ® 1c where U acts on 
the system and bath. Note that since U conserves energy, U is block diagonal in the energy eigenbasis belonging to 
different energies. More precisely, if the total Hamiltonian H = Hs+Hb is block diagonal H = E n T\.E n where n 
is the projector onto the subspace of energy E n = nhco spanned by |0) s |0) B for n = 0 and {|0) s | n) B , |l) s |n — 1) B } 
for n = 1,2,3,..., then U = ©„ Ue „ where Ue„ is a unitary acting only on the subspace of energy E n . That is, 
n Bii t/ E n ^E n = Ue „• 

Consider the unitary transformations Ue.„ defined by the following action: 

Ue 0 |0} s |0 ) b = |0) s |0 ) B =: |^ 0 > , (45) 

u E n |0) s | n) B = Vb\0) s \n) B + y/l -b\l) s \n-l) B =: |^ Bn ) for n = 1,2,3,... , (46) 

UE n |l) s \n-l) B = Vl~b\0) s \n) B - \/fe|l) s \n- 1> B =: |^ Bn > for n = 1,2,3,... , (47) 

where 0 < b < 1 is a parameter that will be chosen in accordance with the desired target state ps below. It will be 
useful to observe that in the subspace {|0) s \n) B , |l) s \n — 1) B }, the unitary Ue „ can be written as 

(48) 

which makes it easy to see that U = W is Hermitian. Note that the states are normalized and ) = 0 for 

n = 1, 2,3,.... The bath thermal state is 


I 00 

TB=^Y, Q ~ nEsP I»X"|b, (49) 

B n—0 

where Zb — Y^=s e ~ nEs ^ = 4/(1 — e~ Es &) is the partition function of the bath, and we have used the fact that 
E n = nhco = nEs . The unitary thus transforms the overall state as 

oo 

U(\0)(0\ s ® t b )U* = — J2 e~ nEsf) U (|0) s | n) B s (0| B (n\) (50) 

B n—0 

1 °° 1 

= -^J2 e ~ nEs0U (\°)s\n) B s (0| B (n|)C/ t + —C/(|0) S |0) B s (0| B (0|)t/t (51) 

n =1 

i 00 i 

= + ^|0)(0| s SS |0)(0| B (52) 

n =1 

=: p°sb • (53) 

By linearity of the partial trace operation, we have that 

TtbG o° sb ) = (6|0)(0| s + (1 - b)\l)(l\ s ) + ^|0)(0| s (54) 

= Po|0)(0| s +p 1 |l)(l| s , (55) 

where 

Po = ^((Z B -l)b+l) , 

Pi = 1 ~ Po ■ 


Note that since 0 < b < 1, po £ [1/Z B , 1] = [1 — e Es P, 1]. Solving (56) for b gives 


(56) 

(57) 








2. 


The reversal operation 


Let us now construct the reversal map lZ a ^ p . Since no catalyst is needed, we can write the reversal map as 

7 Za^pips) = Tr B [UHPS ® t b )U] = Tr B [U(p s ® t b )U*] , (59) 

where we have used the fact that U = U'. To evaluate the reversal map for arbitrary pg, let us first note that by a 
calculation similar to the above 

oo 

[/(|l)(l| s ®f B )[/t = — Ee“ nB ^l^„ +1 X^k +1 l =: Psb ■ (60) 

B n =0 

Using the linearity of the partial trace, we furthermore observe that 

^b[p\b] =(1-6)|0)<0| s + 6|1X1|s - (61) 

Using (53) and (60) together with (55) and (61), we then have 

TZ a ^ p {ps) = Tr s [U(p s <S t b )U^\ (62) 

= PoTrs [p SB ] +PiTr s [p SB ] (63) 

= Po (po|0)(0|s +pi|l)(l|s) +pi ((1 - &)|0)(0|s + 6|l)(l|s) (64) 

= p?\ms+p?m\s, ( 65 ) 

with 

P? ■= 1 -P?, ( 66 ) 

P^ ■= {Pof + (Pi) 2 7 Zb t = (Po) 2 + (1 - Po) 2 e Bs/3 , (67) 

where we have used the fact that Po +Pi = 1 and Z B = 1/(1 — e~ Es ^). We can now compute the lower bound for 
Wi"t no (p -> o). We find 


W?™°(p ^ a) > -kT\ogF{as,n^ P (ps)) (68) 

= -kT log . (69) 

Plugging in (66) followed by (58) into (68) we find 

W.™°(P^ a) >-kT log [(p 0 ) 2 + (l-p 0 ) 2 e^] , (70) 

where we recall p 0 £ [1 /Z B , 1] = [1 — e~ Es ^, 1]. 


3. Three special cases 


We examine three special cases of (70): 

1) Consider p 0 = 1. In this case we want to form the state |0)(0|g from the state |0)(0|s. The work invested must 
clearly be zero in this case. The RHS of (70) is also zero, and hence the bound (17) is tight for this case. 

2) Consider p 0 = 1/Zs = 1/(1 + e~ Es/3 ). That is, we want to form a pure state from the thermal state ps = ts- 
In this case, the RHS of (70) simplifies to — kT log F(as,Pa^ P {ps)) = kT log Zg- By direct calculation using 
the 2nd laws (using Eqs. (13)-(14)) we find lUga// 0 = W-™ no = {log Zg)/ft and thus the bound is also tight for 
this case. 

3) Consider po = 1/Z B . That is, we want to form a pure state from the state whose ground state population 
is the same as the ground state population of the harmonic oscillator bath. In this case, (70) reduces to 
—kTlog F(a s , lZ a ^p(ps)) = -kT log[l + e~ 2E - e~ E ^}. 
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FIG. 1. We show how the different classes of resource theoretic thermodynamic operations are included within each other. The 
most general ones are the maps that preserve the Gibbs state of a system, known as Gibbs Preserving maps (GP), followed 
by those in which we allow for an arbitrary number of catalysts (CCTO, and the slightly wider darker region where we allow 
the catalysts to correlate with the system, nCTO), followed by catalytic thermal operations (CTO) with a single catalyst, and 
finally the simplest one, thermal operations (TO). If one has access to the full set of CCTO, then the standard second law 
determines when transitions are possible, and hence Eq. (8) holds. In the main section, we show that the recovery map of an 
operation within one of the classes is inside the same class, and thus if one can implement all the operations of a specific class, 
they can also implement the corresponding recovery operation without the need for additional resources 


V. EXTENDING TO MORE GENERAL OPERATIONS INVOLVING CATALYSTS 

In this section, we show how the results of Section III hold for increasingly general sets of resource theoretic 
thermodynamic operations. In particular, we give proofs for classes of thermodynamic operations consisting of: 

1. Catalytic thermal operations (CTO), with a single catalyst in a state that is left unchanged at the end, as 
introduced in [20]. 

2. The regime of approximate catalytic thermal operations with various catalysts, where their local states do 
not change, but can become correlated [10]. We refer to this class as “correlating catalytic thermal opera¬ 
tions” (CCTO). 

3. Operations with various catalysts where they also become correlated with the system, which is a potentially 
physically relevant extension that has not been previously defined. We will call this more general class n-catalytic 
thermal operations (nCTO). 

4. The most general kind of maps that can be taken as free operations of a thermodynamic resource theory [42]: 
the set of maps that leave the thermal state invariant, called Gibbs-preserving maps (GP). A resource theory 
defined by these free operations fits the requirements of [43]. 

We show how each set is included in the others in the diagram of Fig 1. 

For 3), it is known [10] that even in the nano regime the free energy is recovered as the unique measure for 
thermodynamical processes. This case is hence of particular importance in this context, as the main results of the 
present work are written in terms of the free energy: if one has access to the full set of such operations, then we are 
in a regime in which Eq. (8) holds. 

First we prove the following lemma, which highlights a particular feature shared by most of maps for which the 
main results can be shown. For the last, more general, set of operations, Gibbs preserving maps, we give a separate 
proof in Section VI, as the result differs in form slightly. 

Lemma 2. Let T(-) be a quantum channel with a full-rank steady state ts = T{ts ) and a Stinespring dilation 

T(-)s = Ti- e [U((-) s ®Pe)U% (71) 

for some unitary U and an environment state p E , such that 

U(t s <S> p E )U t =t s ® p E . 


(72) 
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That is, at the fixed point no correlations with the environment are created. It then holds that for an arbitrary initial 
state ps, and as = T(ps) 

D(ps\\rs) - ^(o-sIIts) > D(p s \\TZ(a s )), (73) 

where IZ(-) is the Petz recovery map, given by 

n-) = Tr E {UH(-) s ®p E )U}. (74) 

Proof. Our proof follows the steps of the particular case of thermal operations shown previously. We first write 

D{ps\\rs) = D(p s 0 p E \\rs 0 pe) (75) 

= D(U(ps®Pe)U^\\U(t s ®pe)U f ) (76) 

= D(U{ps 0 Pe)U^\\ts 0 Pe), (77) 

where we have used the main assumption of the lemma in (72) and the facts that the relative entropy is invariant 
with respect to tensoring an ancilla state or applying a unitary. 

Now we recall the identity of (26) from the proof of Theorem 1: 

D{vcd\\0cd) - D(r] D \\d D ) = Tr(?? cr ,[logr?cD - log 9 C d - log I c 0 Vd + logic ® 0 D }) , (78) 

where suppler?) C supp(0cr>)- We use it together with (75) to write 

D(ps\\rs) - £>(cts||ts) = Tr(17 (ps 0 p E ) U ] [log U(p s 0 p E )U ] - logrg 0 p E ~ log0 Ie + logT S 0 Ie])- (79) 

The last three terms on the right-hand side above can be simplified significantly 

- log t s 0 Pe log a s 0 I E + log t s 0 I E = - log I s 0 p E - log a s 0 I E (80) 

= -logo’s ®PE, (81) 

which leads to 

D{ps\\rs) - D(a s \\T S ) = D(U(p s 0 p E )U^\\a s 0 p E ) (82) 

We also have that 

D{U{ps 0 p E )U ] \\a s 0 Pe) = D(ps 0 p E \\U\as 0 Pe)U). (83) 

Putting everything together, we see that 

D(ps\\ts) - -D(<js||ts) = D(p s 0 p E \\U ] {a s 0 pe)U) (84) 

> D(ps\\nas)). (85) 


What is left is to show that 7 Z is indeed the Petz recovery map. Again this follows through the same reasoning as it 
does for thermal operations. 

The adjoint of the map T(-) is as follows 


(0s-> Trs [p¥uH(-) s ®Ie)Up)' 2 ' 
by definition the Petz recovery channel is given as 


a 


1/2 




Pe 2 U + (r s 1/2 (-)s T s 1/2 ® i e) Up 


1/2 

E 


M 2 


By a series of steps similar to those shown previously, we have that 


i / 2 


Tr £ [a /2 t/l (r" 1/2 (-) s rs 1/2 0 I E ) Up\ 


/2 


1/2 


= Tr K 


(t S ® Pe ) 1/2 U ] (t s 1/2 (-) s t s 1/2 0 I E ) U (t s 0 Pe ) 1/2 
= Tr ® C/l (t s 0 Pe ) 1 ' 2 (ts 1/2 (■) s t s 1/2 ® (ts 0 p ^) 1/2 U 
= Tr E [UH(-) s ®Pe)U] 

= ^(')» 


( 86 ) 

(87) 


( 88 ) 

(89) 

(90) 

(91) 


We note that the matrix (ts 0 pe) 1 ^ 2 is diagonal in the same basis as ts 0 pe, and both have the same degeneracies. 
Hence, if the latter commutes with U, so does the former. This allows us to go from the second to the third line by 
using the assumption of (72). □ 
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This lemma implies that for any quantum channel that has a dilation satisfying the condition in (72), we can show 
a generalization of (34). In the next lemma, we develop a further relevant set of maps for which the condition in (72) 
holds. We call these n-catalytic thermal operations (nCTO), which are defined through an energy-conserving unitary 
V acting on the system S , a set C of n catalysts C = ®” =1 C,; with initial states 77 c., and the bath B, such that 

Tr s [V (ps ® t b ®" = 1 rj Ci ) V*] = &SC , (92) 

where Tr c [asc] = erg, and we impose the constraint T^sc\Ci [^sc] = he, V i, meaning that the local states of the 
catalysts return to their original form. The unitary conserves the energy of bath, system, and all the catalysts, so 
that [V, Hg + H b + Y^i=i He, ] = 0, but correlations between the different parts may be allowed in the final state. We 
can then show the following: 

Lemma 3. For a given n-catalytic thermal operation (nCTO) of the form in (92), if the input state ps is the thermal 
state ts (the fixed point), then 

V (ts 0 T B ®" = 1 77 a) = ts ® t b ®" =1 pa ■ (93) 

Proof. Let p B = Trgqffo (ts ® t b ®" =1 pa) V"l] denote the local state of the bath after the transformation, and denote 
the total Hamiltonian as H = Hg + H B + ^Ci , the sum of all the local ones. Conservation of energy before and 
after the operation corresponds to the following: 

n 

Tr [HV (t s ® t b <g >" =1 r)d) V' 1 } = Tr[iJ(r s ® t b ®>?=i VcJ] = Tr[77 s r s ] + Tt[H b t b ] + ^ Tr [H Ci pCi]- (94) 

2=1 

Also, the total average energy is the sum of the local energies 

n 

Tr [HV (t s 0 t b ®? = i r icJV*] = Tr [H s t s } + Tr [H B p B ] + ^ Tr\H Cl Vc,}, (95) 

and hence Ti \H B p B ] = Tr [H b t b \. For a given fixed energy, the thermal state is the unique state with the highest 
entropy [44, Theorem 1.3], and hence S(t b ) > S(p B ). 

Now we consider the entropy S(p ) = —Tr[plogp] before and after the transformation. Since the joint operation is 
a unitary, we have from unitary invariance and subadditivity of quantum entropy that 

n 

S(t s ) + S(t b ) + ^2 s ( r lGi) = S(v (ts ® T B ®f =1 Pa) y ] ) (96) 

2=1 

n 

< S(t s ) + S(p B )+ ^2S(pc z ), (97) 

i = l 

and hence that S(f B ) < S(p B ). Given our conclusion in (94) regarding conservation of energy, it must then be the 
case that S(t b ) = S(p B ), and thus that the inequality in (97) is an equality. Saturation of subadditivity of entropy 
in (97) occurs if and only if 

V (t s ® t b ®" =1 pa ) V* = t s ® tb ®Li VCi , (98) 

completing the proof. □ 

Putting together Lemmas 2 and 3 and taking the environment state pe from Lemma 2 to be the state of the bath 
and the set of catalysts (i.e., p B = t b ®” =1 pcf), we arrive at the following conclusion: 

Theorem 2. Let T be an nCTO of the form in (92) and given by 

T(-) = Tr BC [U((-)s®f B ®? =1 Pa)U ] ]. (99) 

Then it obeys the inequality 

D(ps\\rs) - £>(o-s||rsr) > D(p s \\'R.(a s )), ( 100 ) 

where IZ(-) is the Petz recovery map, given by 

U(-) = Tr sc [[/t((.)s ® TB ®? = 1 VcM (101) 

which is a thermodynamic operation that belongs to the same class as the original operation T (so for instance, the 
recovery map of a CTO is also a CTO), but with the global reversible dynamics generated by inverse unitary U' 
instead. 
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There are three relevant subsets of ?r-catalytic thermal operations to which Theorem 2 applies: 

1. Catalytic thermal operations, as defined in [20], where we have a single catalyst that does not become correlated 
with the system. 

2. Catalytic thermal operations with multiple catalysts that may become correlated between themselves, as defined 
in [10]. 

3. Operations where the catalysts and the system are allowed to correlate. This set, although slightly more general 
and possibly powerful, has not been studied in detail up to date. 

As already pointed above, the second kind is of particular interest because it can be shown that if one uses 
three or more catalysts in a process in such a way that they correlate between themselves, while their local states 
remain unchanged, the criteria for when thermodynamical transitions are possible reduces to the standard free energy 
decreasing [10], as in (5). 

We have shown that in this context, the Petz recovery map can be understood as the map generated by the inverse 
of a global dynamics, which can be understood as a time reversal of the original map. This fact has been noted in 
quantum thermodynamics in the context of fluctuation theorems [45-47]. 


VI. GIBBS PRESERVING MAPS 

The only set of maps to which the conditions of Lemma 2 do not apply is the most general set of Gibbs preserving 
maps [42], and we hence need a different method. To prove a bound of the form of (21), we use the following general 
result for quantum maps: 

Theorem 3 (From [48]). Let Nf) be a quantum channel, and let r) and 6 be quantum states. We have that 

D(vP) - D(Af( V )\\Af(9)) > -2log F(p, [ dt p(t) Af t (Af(y))), (102) 

Jr 

where F(p,a) = Tr[\/ \fop\fa\ is the quantum fidelity, the map Aft is the following rotated recovery map 

Aft(-) = e lt/ 2 Af{Af{e)- it / 2 {-)N{e) lt/ 2 )e- it '‘ 2 , ( 103 ) 

with Af defined as in (38) and p(t) = ^(cosh(7r t) + l) -1 is a probability density function. 

It can be seen by inspection that if we take the map Af to be Gibbs-preserving so that ts = Af(rs) and if we set 
9 = ts, then the rotated recovery map is Gibbs-preserving as well, namely Aft(rs) = ts- More explicitly, the bound 
on A is as follows: 

A = D(ps\\t s ) - T>(<ts||ts) > -21ogF(p s , [ dt p(t)Af t (a s )), (104) 

Jr 

where instead of having the relative entropy, we have the fidelity in the lower bound for the decrease of free energy. 


VII. CONCLUSION 

We have shown that the amount of work drawn from a process is directly linked to the reversibility of the process. 
Specifically, we see that there are some regimes where if the amount of work drawn is small, then there exists a 
recovery operation that approximately restores the system to its initial state at no work cost at all. It is interesting 
to note that our result also makes a statement about the standard second law whenever ps can be transformed to as- 
When A is small, then the transformation ps to as can be reversed very well by the map TZ p ^ a . 

Our main result applies to all regimes where the standard free energy is relevant, and it is a very interesting open 
question to extend our result to regimes in which we require the full set of second laws [20]. What makes this question 
challenging is that (26) does not carry over to the regime of D a for a ^ 1, and indeed recent work [49] suggests 
that other quantities naturally generalize the difference of relative entropies—and this generalization does not always 
result in the difference of a-Renyi relative entropies. It hence forms a more fundamental challenge to understand 
whether the difference of such a-relative entropies, or the quantities suggested in [49] should be our starting point. 
However, the quantities in [49] would require a proof of a new set of second laws. 
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We have applied our analysis to the case of investing work, which in the regime where only the standard free energy 
is relevant can be characterized fully by how much work can be gained by the inverse process. This relation to the 
inverse process is not true in the nano-regime where all the refined second laws of [20] become relevant. Nevertheless, 
we have shown that the reversal operation of said inverse process can indeed be used to understand the amount of 
work that needs to be invested, adding another piece to the growing puzzle that is quantum thermodynamics. 

Since this initial work, there has been a series of recent results consisting in giving lower bounds to the decrease 
of relative entropy in different cases of interest, covering a number of different branches within quantum information 
theory, such as [50-54], which shows the importance of the concept of recoverability and of recovery maps. 
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